arXiv:1502.01336v2 [gr-qc] 5 May 2015 


Renormalized vacuum polarization of rotating black holes* 

Hugo R. C. Ferreira^ 


School of Mathematical Sciences, University of Nottingham, 
Nottingham NG7 2RD, United Kingdom 
pmxhrf Qnottingham.ac.uk 


Abstract 

Quantum field theory on rotating black hole spacetimes is plagued with technical dif¬ 
ficulties. Here, we describe a general method to renormalize and compute the vacuum 
polarization of a quantum field in the Hartle-Hawking state on rotating black holes. We 
exemplify the technique with a massive scalar field on the warped AdSa black hole solu¬ 
tion to topologically massive gravity, a deformation of (2-|-l)-dimensional Einstein gravity. 
We use a “quasi-Euclidean” technique, which generalizes the Euclidean techniques used for 
static spacetimes, and we subtract the divergences by matching to a sum over mode solu¬ 
tions on Minkowski spacetime. This allows us, for the first time, to have a general method 
to compute the renormalized vacuum polarization (and, more importantly, the renormal¬ 
ized stress-energy tensor), for a given quantum state, on a rotating black hole, such as the 
physically relevant case of the Kerr black hole in four dimensions. 


1 Introduction 

The study of quantum field theory on stationary black hole spacetimes has proven to be much 
more challenging than on static black hole spacetimes. In particular, it has not been possible 
to compute the expectation values of the renormalized stress-energy tensor for a matter field 
in a given quantum state on the Kerr spacetime. Apart from the technical complexities of the 
calculations, other difficulties include the nonexistence of generalizations of the Hartle-Hawking 
state and of Euclidean methods used in static spacetimes. 

Recently, Ref. laid down a method to compute the renormalized vacuum polarization of 
a quantum field in the Hartle-Hawking state on a rotating warped AdSs black hole surrounded 
by a mirror with Dirichlet boundary conditions, implementing the Hadamard renormalization 
prescription on the complex Riemannian section of the spacetime. Here, I review this method 
and extend the results for a generic (2-|-l)-dimensional rotating black hole spacetime. This 
calculation can be taken as a warm-up for the computation of the more physically relevant 
renormalized stress-energy tensor. We anticipate that this technique can be extended to a wider 
class of rotating black hole spacetimes, in particular the Kerr black hole in four dimensions. 

The contents of this note are as follows. In Sec. [2]|^ we outline the method for a generic 
(2-|-I)-dimensional stationary black hole spacetime and apply it in Sec. to the case of the 
warped AdSs black hole. We present our conclusions in Sec. and the proof of a theorem is 
left to Appendix]^ Throughout this note we use the signature and units in which 

h = c = G = kB = l- 

‘Based on a talk given by H.R.C.F. at “VII Black Holes Workshop”, Aveiro, 18-19 December 2014. 

^Work done in collaboration with Jorma Louko. 
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2 Scalar field on a rotating black hole spacetime 


For concreteness, in the following we will consider a generic (2+l)-dimensional stationary black 
hole spacetime, whose metric is of the form described below. This is the case of the warped 
AdSs black hole described in Sec. However, as argued, the method should be applicable to 
a wide range of rotating black hole spacetimes in three and more spacetime dimensions. 

In spherical coordinates the metric of a (2+l)-dimensional stationary black hole 

can be written as 

ds'^ =-N^{r) dt'^ + grrir) dr"^ + g 0 g{r) (^d9 + N^{r) dt^ , (1) 

where N{r) is the lapse function and N^{r) is the shift function. Let the event horizon be at 
r = r_|_ and let be the angular velocity of the horizon with respect to the coordinate system. 
Denote by x = dg the Killing vector field which generates the horizon and which is 

timelike for r G (r_|_,rc), where r = rc is the location of the speed-of-light surface, if it exists. 

For most cases of physical interest, such as the Kerr spacetime, there is not any timelike 
Killing vector field in all of the exterior region r > r+. As a consequence, in the context of 
quantum field theory, there is not a well defined quantum vacuum state which is regular at the 
horizon and is invariant under the isometries of the spacetime [^. However, a vacuum state 
with these properties can be defined if we restrict the spacetime by inserting a boundary Ai 
at constant radius r = r_M, with r_|_ < rjn < rc, in which the field satisfies Dirichlet boundary 
conditions, as in this region x is a timelike Killing vector field. 

For convenience, we change to “corotating coordinates” {t = t, r, 6 = 6 — such that 

X is given hy x = and the metric is then given by 

ds'^ =-N‘^{r) dP + grrir) dr^ + ggo{r) (^dO + {N^{r) + Q-}i)dtj . (2) 

From now on, we only consider the portion of the exterior region from the horizon up to the 
mirror, which we denote by I, in contrast with the full exterior region I. A real massive scalar 
field <I> on I obeys the Klein-Gordon equation 

(V^ -m'^ -CR)^ = 0, (3) 

where m is the mass of the field, R is the Ricci scalar and ^ is the curvature coupling parameter. 
We consider mode solutions of Q of the form 


where a) G M and A; G Z. 

Given the construction above, it is possible to show that there is a natural positive and 
negative frequency decomposition of the mode solutions with respect to the affine parameters 
on the horizon [^. It is then possible to canonically quantize the scalar field and define a 
vacuum state \H), which is regular at the horizons and invariant under the spacetime isometries. 
Therefore, \H) has all of the properties of the “Hartle-Hawking vacuum state”, as usually defined 
for static black hole spacetimes [^. 

The Feynman propagator is defined as 

G^{x,x') := i {H\^ {^x)<^{x')) \H), (5) 
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where is the time-ordering operator and <h(x) is interpreted as an operator-valued distribution 
which acts on the Fock space of the theory. The Feynman propagator is one of the Green’s 
functions associated with the Klein-Gordon equation. 

The Feynman propagator, evaluated for certain quantum states states, has a Hadamard 
expansion of the form 


{x, x') = 


[/(x, x') 


4\/27r I ^Jcr{x, x') + ie 


P -|- W{x,x') 


0 + . 


( 6 ) 


Here, x and x' belong to a neighborhood N and are linked by a unique geodesic which lies 
entirely in N-, a{x,x') is half of the squared geodesic distance between x and x'; U{x,x') and 
W{x,x') are symmetric and regular biscalar functions. Note that this expansion is only valid 
for (2-|-l)-dimensional spacetimes. Spacetimes with d > 2 dimensions have a singular term 
of the form [a{x,x') + and, in the case of spacetimes with even dimension d > 4, an 

additional singular contribution of the form log[cr(x, x') -|- ie\ is present. 

A quantum state for which the short-distance singularity structure of is given by Q (or 
its higher-dimensional versions) is called a “Hadamard state”. Examples include the Minkowski 
state on the Minkowski spacetime and the Hartle-Hawking state on the Schwarszchild black 
hole [^. Even though, as noted above, there is not an analogue of the Hartle-Hawking state 
on a rotating black hole [^, an isometry-invariant state which is regular at the horizons can be 
defined if a Dirichlet boundary is present, so that the quantum field cannot reach the speed-of- 
light surface. We still call this state the “Hartle-Hawking state”. 

Goncerning the regular functions U{x,x') and VF(x,x') in the expansion ([^, it can be 
shown (see e.g. Ref. 0) that U (x, x') only depends on the geometry along the geodesics joining 
X to x', whereas W{x,x') contains the quantum state dependence of the Feynman propagator. 
Therefore, the singular, state-independent part of the Feynman propagator, the “Hadamard 
singular part”, is 


G'Had(a:,x') 


i U{x,x') 
4\/27r ^Ja{x] x') + ie 


( 7 ) 


The renormalized vacuum polarization (<h^(x)) in any Hadamard state is then defined to be 


(^2(x)) .= lim \p^[x,x) - GHad(a:,x')] • 
x'^x 


( 8 ) 


At this stage, the Euclidean techniques used for static spacetimes to more easily compute 
the Feynman propagator are not available for stationary, but not static, spacetimes. In the 
next section, we describe a “quasi-Euclidean” method to achieve the same result. 


3 Quasi-Euclidean method 

Euclidean methods can be a powerful tool to do calculations in quantum field theory on static 
spacetimes. A static spacetime can be thought of as a real Lorentzian section of a complex 
manifold, for which it is always possible to find a real Riemannian (or “Euclidean”) section by 
performing an appropriate analytical continuation. For a (2-|-l)-dimensional static spacetime 
whose metric in coordinates (t, r, 9) is 

ds^ = —N‘^{r) dt^ + grr{r) dr^ -|- geeir) dO^ , (9) 
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where t is a global timelike coordinate, one can obtain the real Riemannian section by performing 
a Wick rotation t —)• —ir, r G M, 

= N‘^{r) + grr{r) dr^ + geeir) dO"^ . ( 10 ) 

This analytic continuation procedure does not easily generalize to spacetimes that are sta¬ 
tionary but not static. For instance, for the exterior of a rotating black hole, there might not 
be a globally defined timelike vector field or an analytic continuation in the coordinates that 
results in a real section with a positive definite metric. Both of these issues are present in 
Kerr [^. 

If we only consider region I of the (2-|-l)-dimensional rotating black hole spacetime, there 
exists an everywhere timelike Killing vector field, x = If we now perform a Wick rotation 
i = —ir, with r G M, the metric Q becomes 

ds'^ = N‘^{r) dr'^ + grrir) dr"^ + gee (^d9 — i [N^{r) + ^l'}i)dT^ . (11) 

This is the complex-valued metric g'^ of the “complex Riemannian” (or “quasi-Euclidean”) 
section of a complex manifold, in which region I is a real Lorentzian section [t-IO . This 
metric is regular at the horizon if r is periodic with period 27r/K_|_, where is the surface 
gravity of the black hole. 

The Green’s function G associated with the Klein-Gordon equation in the complex Rieman¬ 
nian section satisfies the distributional equation 

6^{x,x') 5 {t — T')5{r — r')5{6 — &) 


(V^ — w?) G{x, x') = -- 






( 12 ) 


where g{x) := \ det( 5 ^;,)| and := 

In contrast to the real Lorentzian section, there is a unique solution to this equation in the 
complex Riemannian section which satisfies the following boundary conditions: (i) G{x, x') is 
regular at r = r+, and (ii) G{x,x') satisfies the Dirichlet boundary conditions at r = This 
is because two of the directions of the complex spacetime are periodic, while the third direction 
is compact. In contrast, on static spacetimes without any boundary (and suitable asymptotic 
properties at infinity), the Euclidean section has a unique Euclidean Green’s function, due to 
the ellipticity of the Klein-Gordon operator. 

Given the periodicity conditions of r and 9, one has 


5{t — t') 

_ ^ iK+n{T-T') 

27r ^ 

n=—oo 

(13) 

5{9 - 9') 

^ oo 

_ I ik{e-0') 

2tt ^ 

(14) 


k=—oo 


understood as distributional identities. We can than find the unique solution of Eq. (12) and 
write it as a mode sum. 




CXD 

E ■ 

n=—oo 


^iK,+ n{T—T') 


^ Gnk{r,r') =:Y,Gnk{x,x') 


(15) 


k=—oo 


n,k 


In practice, the radial part of the Green’s function GnkY, r') is found by solving the radial 
field equation, using standard Green’s function techniques. In 2-|-l dimensions it is generally 
possible to write it in closed form, whereas we need to resort to numerical methods for four or 
more dimensions. 
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4 Hadamard renormalization 


As with with the Feynman propagator in the real Lorentz section, the Green’s function G in 
the complex Riemannian section has an Hadamard expansion and, in particular, we can write 
its Hadamard singular part as 


G^Had(^) ^ ) 


1 U{x,x') 

4-v/27r a{x, x') ’ 


(16) 


where (t(x, x') is the generalization of the quantity defined in the real Lorentz section (see 
Ref. for more details). 

In an analogous way to the Lorentzian case, we now subtract the Hadamard singular part 
from the Green’s function G and obtain the vacuum polarization. 


^ [G(x,x') - GHad(a:,x')] • 


(17) 


By construction, the Green’s function G is regular at r = r_|_, satisfies the Dirichlet boundary 
conditions at r = r _\4 and is invariant under the spacetime isometries. Therefore, (<l>^(x)) as 
given by 0 is the vacuum polarization for a scalar field in the Hartle-Hawking state. 

It remains to perform the subtraction in ( |17[ ) before the coincidence limit can be taken. As 
G is known only as the mode sum (15), the evaluation of (<I>^(x)) requires Guad to be rewritten 
as a mode sum that can be combined with ( |15[ ) so that the divergences in the coincidence 
limit get subtracted under the sum term by term. We accomplish this by comparing Guad to 
the Hadamard singular part for a scalar field in rotating Minkowski spacetime in the complex 
Riemannian section, and by rewriting it as a sum over Minkowski mode solutions. 

It can be shown that the Hadamard singular part of the Green’s function for a scalar field 
in the Minkowski vacuum can be written as 


n,k 


nk 


(x,x') + G^g(x,x') 


(18) 


where the first term on the rhs is the full Green’s function and Gjfg(x,x') is a term which is 
finite when x' —)> x [^. 

It is convenient at this stage to consider a particular choice of point separation. Assume 
that the black hole metric is given in coordinates (r, r, 0), whereas the Minkowski metric is 
given in coordinates (r, p, 0). Now, consider the case of angular separation in each spacetime, 
such that for the black hole case x = (r, r, 0) and x' = {T,r,6), with 0 > 0, and similarly for 
the Minkowski case. 

The expansion of the Hadamard singular parts for small 9 are 


GK>^d{x,x') = 
G^^{x,x') = 


^ 1 + 00), 


Vmi0 ^ 


dvr 




^ + O0), 
9 


(19) 

( 20 ) 


where 5'||(/o) = is the 00-component of the metric for the rotating Minkowski spacetime. We 
are free to identify 


g^0p)^l ‘^{r)gee{r), 


( 21 ) 
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where 7 (r) > 0 is a function to be specified. This identification provides a matching between 
the two radial coordinates, p = p{r) = 7 ~^(r)Y^ 5 gg(r). 

Given this identification, we can now write 


G[x,x') - GHad(a:,x') = ^ \G^f{x,x')--f ^G^,,{x,x') 

n,k 


.-1 


■7 


(x,x') + O(0). 


( 22 ) 


The Minkowski Green’s function has several free parameters; p (radial coordinate), Tm 
(temperature of the scalar field), (angular velocity of the coordinate system) and 
(squared mass of the scalar field), besides the unspecified factor 7 we introduced above. These 
can be chosen such that the double sum in (22) is convergent when 0 — 7 - 0 . 


The main result of this section is the following: 

Theorem 4.1. If the parameters 7 , Tm and Om are chosen as 


ry[r) = N{r), Tm = — , Vim = {r) + Hh ■, 

zvr 


then the double sum in 


Proof. See Appendix 


is finite in the coincidence limit. 


(23) 


□ 


This choice corresponds to have the temperature Tm of the scalar field in Minkowski to 
match the Hawking temperature of the black hole and to have the angular velocity to 
be equal to the one measured by a locally nonrotating observer at radius r in the black hole 
spacetime. 

The key aspect of the proof of this theorem is that, in order to remove the divergences, we 
only need to know the asymptotic behaviour of the Green’s function summands G™(r, r) and 
G^(p, p) for large values of the quantum number n and k, and not the full solutions. This 
implies that, apart from technical difficulties, this method can be applied to black holes in four or 
more dimensions, for which although we can only obtain the Green’s functions numerically, the 
asymptotic expansions of the summands for large quantum numbers can be explicitly computed 
using the procedure described in Appendix [A| (and in Ref. in more detail). 

Setting the parameters as in ( pS] ), it is now possible to take the coincidence limit 0 — 7-0 
of (22) and compute the renormalized vacuum polarization ( |17[ ). In the next section, as an 
example, we present the numerical results for the particular case of the warped AdSa black 
hole. 


5 Numerical results for the warped AdSa black hole 

In this section, we exemplify the method described above by numerically computing the vacuum 
polarization of the scalar field in the Hartle-Hawking state for a (2+l)-dimensional rotating 
black hole, the spacelike stretched black hole. This is one of the several types of warped AdSa 

to topologically massive gravity , an extension of Einstein gravity 


black hole solutions 11 


in 2+1 dimensions with a propagating degree of freedom. Its metric, in coordinates (t, r, 0), is 
given b 3{3 


= -N‘^{r)dt^ + 


dr"^ 


4R2(r-)Ar2(r; 


+ R^{r) (^d9 + N^{r 


(24) 


^We set the cosmological length £ = 1 . 
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Figure 1: Carter-Penrose diagrams of the spacelike stretched black hole spacetime for the case 
ro < r_ < r_|_ on the left and a portion of the spacetime with the mirrors in place on the right. 


with 


R^{r) 

N'^{r) 

N^{r) 


- - l)r + + 3)(r+ + r_) - 4j/v^r+r_(z/2 + 3) , 

+ 3)(r — r_|_)(r — r_) 

4i?2(r) ’ 

2 i/r — -yr_|_r_(z^2 _|_ 3) 

2i?2(r-) ■ 


(25a) 

(25b) 

(25c) 


There are outer and inner horizons at r = r+ and r = r_, respectively, where the coor¬ 
dinates {t,r,6) become singular, and a singularity at r = ro (the largest zero of i?^(r)). The 
dimensionless coupling v G (l,oo) is the warp factor, and in the limit u ^ 1 the above metric 
reduces to the metric of the BTZ black hole 


a speed-of-light surface is located at 


14,15 in a rotating frame. For > 1, there exists 


r = rc = 


4z^^r+ — (z/^ -|- 3)r_ 
_ 


(26) 


We assume that a mirror A1 is present at constant radius r = r_M < ^C) at which Dirichlet 
boundary conditions are imposed. Moreover, the angular velocity of the outer horizon is given 
by 


= 


2 

2z/r+ — i/r+r_(z/2 -|- 3) 


(27) 


Note in particular that this spacetime does not have a static limit. 

The Carter-Penrose diagram for this spacetime when rg < r_ < r+ is shown in Fig. which 
is essentially of the same form of those of asymptotically flat black holes in 3-I-1 dimensions. 
More details about this black hole solution can be found in Ref. 


16 and references therein. 


7 























As described in Sec.j^ we restrict our attention only to region I of the spacetime (see Fig.[y . 
In this case, the vacuum polarization of the scalar field in the Hartle-Hawking state is given by 


OO OO 




k=-oo n=—oo 


yNlr)^ N{r) 




+ 


1 


47rN{r) 


-mm + ^ 


—rriM 


( ^) +*" sgn(C!MiV) 

\/(€) - 1^sinh2(§^) +zesgn(nMiV) 


(28) 


with e —7- 0+ indicating the choice of branch of the square root . 

The numerical results for selected values of the parameters are presented in Fig. In 
the plot, (<h^(x)) is shown as a function of the normalized radial coordinate zjzM^ where 
z = {r — r+)/(r — r_). The plot is very similar to the one obtained in Ref. |17| for a scalar 
field in the (3+l)-dimensional Minkowski spacetime surrounded by a mirror with Dirichlet 
boundary conditions (note that “rotating Minkowski spacetime” is related to “static Minkowski 
spacetime” by a coordinate transformation, hence the results for (<h^(x)) are the same for both 
cases). 

Furthermore, note that (<h^(x)) gets arbitrarily large and negative as the mirror is ap¬ 
proached. This is to be expected, as we imposed that the Green’s function G(x, x') must vanish, 
even when x' —)• x, whereas the subtraction term still diverges when x' —)• x (see chapter 4.3 
of for more details). 



zjzM 

Figure 2: Vacuum polarization for the scalar field as a function of the normalized radial 
coordinate zjzM^ for v = 1.2, r+ = 15, r_ = 1, rj^ = 62 and m = 1. 


6 Conclusions 

In this note, we have described a method to compute the vacuum polarization for a quantum 
field in the Hartle-Hawking state on a general rotating black hole surrounded by a mirror. 



















We have employed a “quasi-Euclidean” technique to obtain the complex Riemannian section 
of the original spacetime, in which the unique Green’s function associated with the Klein- 
Gordon equation can be found. This Green’s function is given as a mode sum and its singular 
behavior in the coincidence limit can be subtracted by a sum over Minkowski modes with the 
same singularity structure. Taking the coincidence limit of this subtraction gives precisely the 
renormalized vacuum polarization. We exemplified the technique with a massive scalar field on 
a ( 2 +l)-dimensional rotating spacelike stretched black hole. 

A key ingredient in our implementation of the Hadamard renormalization was to match the 
mode sum for the Green’s function in the complex Riemannian section of the black hole to a 
mode sum in the complex Riemannian section of a rotating Minkowski spacetime. We anticipate 
that, despite the technical complexity of the calculations, this method can be applied to the Kerr 
spacetime in four dimensions. In this case, the relevant mode solutions to the Klein-Gordon 
equation on the complex Riemannian section would need to be obtained numerically, but the 
asymptotic properties of the solutions for large values of the quantum numbers should be within 
analytic reach, and it is only these asymptotic properties that are required in the matching to 
the Minkowski mode solutions. It should prove interesting to attempt the implementation of 
our method for Kerr in practice. 
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A Proof of Theorem 14.11 

The proof below is a generalization of the proof given in Ref. for the case of the warped 
AdSs black hole discussed in Sec. 

In order to check the convergence of the double sum (22) in the coincidence limit, we obtain 


the asymptotic behaviour of the summand for large values of the quantum numbers n and k. 

For a black hole spacetime with metric ( |11[ ) in the complex Riemannian section, the Klein- 
Gordon equation 

(V2-m2)$(T,r,0) = 0, (29) 

together with the ansatz r, 9) = <^^^( 7 -)^ Rad to 


1 d [ [n+n + ik[N^ + Vty)) k"^ 2 

yisW*' Tr) ---™ 


iV2 


4^nk 0 ' 


Define a new radial coordinate ^ such that the (30) can be written in the form 




(30) 


(31) 


where XnkiO contains all the n and k dependence and is large whenever is large. From 

(30) we obtain 

d L — y'T' 

= * 


9 


(32) 











and 


xlk = 9ee{^^+n + ik{N^+ 0.^)^ + N'^k'^ , rf = g-^^N' 


2 2 
m 


(33) 


Using the results of Appendix D of Ref. the summand G^^{x, x) of the Green’s function 
G{x,x') in (15) has the following asymptotic expansion for large Xnk 




(34) 


Analogously, for the Minkowski Green’s function, the summand G^(x, x) in (18) has the asymp¬ 
totic expansion 


Gfki^,x') = 


1 


2vr 2x: 


+ 0 


nk 


-3 


where 


{p) — (27rrMn. + z/cHm) + k^ ■ 


(35) 


(36) 


The double sum in (22) will be finite in the coincidence limit if the leading term in the 
asymptotic expansion of the summand vanishes, that is, if the term of order x^k expansion 


of G^^{x, x) cancels with the term of order (x^^) ^ of the expansion of 7 ^(r) G^(x, x). This 
only occurs if the free parameters 7 , Tm and Om are chosen as 


j{r) = N{r), Tm = ^ , Qm = A^®(?’) + • 

Ztt 


(37) 


To show that the double sum is indeed finite in the coincidence limit, we need to check that 
the double sum of the remaining terms in the asymptotic expansion of the summand, which 
areO(x;|), is finite. This can be done in the exact same way as in the proof of Theorem III.l 
of Ref. . 
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